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Abstract 



In a recent analysis of the world data on polarized DIS, Bliimlein 
and Bottcher conclude that there is no evidence for higher twist contri- 
butions, in contrast to the claim of the LSS group, who find evidence 
C^ ' for significant higher twist effects. We explain the origin of the appar- 

ent contradiction between these results. 



PACS numbers: 13.60.Hb, 12.38.-t, 14.20.Dh 



1 Introduction 

Because of the limited range of Q^ available in the world data sample for 
polarized DIS, and because of the great accuracy of some of the data at 
relatively low values of Q^, it is important to be able to include the latter data 
in QCD analyses aimed at extracting information on the polarized parton 
densities. To do this consistently the theoretical formulae must be extended 
beyond leading twist (LT) to allow for higher twist (HT) contributions. This 
has been done for some years by Leader, Sidorov and Stamenov (LSS) [1], 
who have shown that such contributions are important, especially in fitting 
the very accurate CLAS data [2]. In contrast, Bliimlein and Bottcher (BB) [3] 
claim that HT effects are completely negligible. In this note we explain that 
the disagreement is probably only apparent, and results from attempting to 
compare two incommensurate quantities. 

2 Structure of higher twist terms 

On the basis of the Operator Product Expansion (OPE), LSS have uti- 
lized an expression for the experimental spin dependent structure function 
gi{x,Q'^)exp of the form 

9l{x,Q'^)exp = 9l{x,Q'^)LT + 9lix,Q'^)TMC + 9lix,Q'^)HT 

= g,(x,Q')LT + 9l{x,Q')TMC+^ + 0{AyQ'), (1) 

where the powers in l/Q"^ corrections to gi have been split into the exactly 
known kinematical target mass corrections (TMC) and genuine dynamical 
higher twist terms (HT). In the LSS analysis the logarithmic Q^ dependence 
of h, which is unknown in QCD, is neglected. Compared to the principal 
l/Q'^ dependence it is expected to be small and the accuracy of the present 
data does not allow its determination. Therefore, the extracted from the 
data values of h{x) correspond to the mean Q^ for each x-bin. The most 
recent results for h{x) for proton and neutron targets [Ij are shown in Fig.l. 
BB write the expression for 9i{x,Q'^)exp in a different form, namely 

9i{x, Q\xp = 9i{x, Q^)lt [I + ^] , (2) 
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Figure 1. Higher twist corrections extracted from the world polarized DIS data 
(corresponding to 5 and 7 x bins). 



where the assumption 



9i{x,Q'^)ht = gi{x,Q^)LT 



C{x) 



(3) 



is used for the twist-4 contribution, and any Q^ dependence in C{x) is ne- 
glected. (It is not clear whether the TMC are accounted for in gi{x, Q'^)lt-) 

BB find no evidence for HT i.e. their C{x) for protons and neutrons is 
compatible with zero. (There is a worrying issue concerning their derivation 
of the neutron value, which we shall comment on later.) 

In trying to understand the apparent discrepancy with LSS, BB write: 
"This result is in disagreement to Ref. [21] (i.e. LSS [5]). Note that in the 
latter analysis a partonic description of Fi(x, Q^) down to low values of Q^ 
is used, while we refer to the measured function." 

This statement originates in the fact that what is actually measured is 
effectively gi{x,Q'^)exp/ Fi{x,Q'^)exp, so that gi{x,Q'^)exp must be extracted 
from the data by multiplying by Fi(x, Q^)exp- 

However, the comment about Ref. [21] (Ref. [5] in this paper) is totally 
incorrect. Indeed in [S] LSS write: 



"According to this method, the gi/Fi and Ai data have been fitted using 
the experimental data for the unpolarized structure function Fi{x,Q'^) 



9i{x,Q" 



91 {x,Q^)lt + 91 {x, Q^)t mc + h{x)/Q^ 



As usual, Fi is replaced by its expression in terms of F2 and R, and phe- 
nomenological parametrizations of the experimental data for F2{x,Q'^) [8] 
and the ratio R{x, Q'^) of the longitudinal to transverse cross-sections [9] are 
used." 

Hence the apparent discrepancy between BB and LSS is certainly not due 
to any difference in the handling of -Fi(a;, Q"^). 

3 Why the discrepancy is only apparent 

On equating Eqs. ([I]) and ([2]) we have 

^^") = a (xt\ ■ (') 

9i[x,(^ )lt 

Thus, if it is legitimate to neglect the scale dependence in h{x) then C{x) 
must vary with Q^, contradicting the use of C{x) as Q^-independent. If, 
on the other hand, it is legitimate to neglect the Q^ dependence in C{x), 
then h{x) must vary with Q^. We thus see that the two approaches are in- 
compatible and their results incommensurate. One of the two methods (or 
perhaps both) has to be incorrect and the fact that their results disagree is 
inevitable and requires no further explanation. However, since the LSS for- 
mulation is closer in structure to the operator product expansion, we believe 
it is more likely to be the correct way to implement HT corrections. More- 
over, the LSS results on HT are in good agreement with those obtained from 
the study of the first moments of the spin structure functions g{''^\x,Q'^), 
and in particular, of the non-singlet structure function gf' (see Ref. [6]). 

One further point remains. BB utilize Eq. (|2]) for proton and deuteron 
data and then extract the neutron value of C{x) via 

Cnix) = ^^Cd{x) - Cp{x) (6) 



where ur, = 0.05 ± 0.01. From Eq. ([2]) one sees that this is incorrect. The 
correct relation should be 



Cn{x) = — -, — TT^T — [- —: — gid{x,Q^)LTCd{x) - gip{x,Q'^)LTCp{x) 

gin\X,Ui )lT J- — i-.OUJD 

(7) 
Thus, even if it is correct to take Cp{x) and Cd{x) independent of Q^, Cn 
will inevitably inherit some dependence on Q^. Note also that the neutron 
spin structure function gin{x,Q'^)LT passes through zero as a function of x 
and it is therefore dangerous to use the above equation to extract the HT 
correction C„. 

4 Conclusions 

We have shown that the LSS and BB methods of extracting HT corrections 
in polarized DIS are incompatible and that it thus makes no sense to compare 
their results - they are incommensurate. We believe that the LSS approach, 
because it is closely related to the operator product expansion, is more likely 
to be the correct one. 
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